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Abstract. Mazur, Tate, and Teitelbaum gave ap-adic analogue of the Birch and Swinnerton-Dyer 
conjecture for elliptic curves. We provide a generalization of their conjecture in the good ordinary 
case to higher dimensional modular abelian varieties over the rationals by constructing the p- 
adic L-function of a modular abelian variety and showing it satisfies the appropriate interpolation 
property. We describe the techniques used to formulate the conjecture and give evidence supporting 
the conjecture in the case when the modular abelian variety is of dimension 2. 

1. Introduction 

The Birch and Swinnerton-Dycr (BSD) conjecture gives a precise relationship between several 
arithmetic invariants of an abelian variety A over a number field K . As formulated by Tate [45], the 
conjecture states the following: 

Conjecture 1.1 (BSD conjecture for abelian varieties). Let A be an abelian variety of dimension g 
over a number field K , and let be its dual. Then the Mordell-Weil rank r of A{K) is equal to the 
analytic rank ovds^i L{A, s) of A and 

lim(. - D-HA, s) . ^A^iA/K)\-Re,iA).U.c 

where Dk is the absolute discriminant of K , U,a is the real period, Reg(^) is the regulator, Cy is the 
Tamagawa number at a finite place v of K , III(^/if) the Shafarevich-Tate group of A and A(A')tors 
is the torsion subgroup of A{K). 

Note that this conjecture relies on two assumptions: that the Shafarevich-Tate group III is finite 
and that the L-series can be analytically continued to s = 1. An analytic continuation is known 
to exist for modular abelian varieties over Q, where an abelian variety is said to be modular if it is 
a quotient of Ji{N) for some level N. In particular, for an elliptic curve E/Q of rank r, the BSD 
conjecture predicts 

Conjecture 1.2 (BSD conjecture for elliptic curves). Let E be an elliptic curve over Q. Then the 
Mordell- Weil rank r of E(Q) is equal to the analytic rank of E and 

Imis-l) L{E,s)^ — — '—. 

\E[Q)torsr 

In 1986, Mazur, Tate, and Teitelbaum [53] gave a p-adic analogue of this conjecture for an elliptic 
curve E over the rationals and a prime p of good ordinary or multiplicative reduction. Much work 
has been done towards a proof of the conjecture, and more is known about the p-adic conjecture 
than its classical counterpart. We give a brief overview of the circle of ideas involved; see also the 
recent work of Stein- Wuthrich [43]. For simplicity, we assume that p > 2. 
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Let ooG denote the Galois group Gal(Q(/ip=o )/Q), where Q(/ipoo) is the cyclotomic extension of Q 
obtained by adjoining all p-power roots of unity, let k : ooG — > denote the cyclotomic character 
and let 7 be a topological generator of F ooG'-^"^-'. For an elliptic curve -E/Q and a prime p 
such that E has good or multiplicative reduction at p, we denote the p-adic regulator, divided by 
logp(K(7))'', by Reg^(i?/Q) and we let Cp{E, T) denote the series expansion of the p-adic L-function 
Lp{E, s) associated to £' in T = K{'yy^^ — 1. 

Conjecture 1.3 (p-adic BSD conjecture for elliptic curves). Let E be an elliptic curve over Q and 
let p be a prime number such that E has good ordinary or multiplicative reduction at p. 

(i) The order of vanishing OTdT{Cp{E,T)) of Cp{E,T) at T ~ is equal to the rank r of E{Q) 
if E has good ordinary or nonspUt multiplicative reduction at p. If E has split multiplicative 
reduction at p, then ovAt{Cp{E, T)) = r + 1. 

(ii) If E has good ordinary or nonsplit multiplicative reduction atp, then the leading term Cp{E,Q) 
satisfies 



(1.1) Cl{E,Q)^ep{E) 



|m(i;/Q)|.Reg^(£;)-n 

|2 



where £p{E) = (l — a^^)'' for a unit root a ofx^ — OpX+p G Qp[x] (with Op the Hecke eigenvalue 
of the newform associated to E) and b is 2 if E has good ordinary reduction at p and 1 if E 
has nonsplit multiplicative reduction at p. 

If E has split multiplicative reduction at p, then 

. ^ |m(i;/Q)|.Reg^(ii;).n.c. 

' ' logp(A^(7)) ■ \E{QU.? 
where ,5^p = o°.j''(^^) o.'nd qE is the Tate period of E over Qp. 

For primes of supcrsingular reduction, the p-adic L-function and the p-adic regulator can also be 
defined, and a p-adic BSD conjecture has been formulated by Bernardi and Perrin-Riou [5]. 

Much work has been done toward a proof of Conjecture 11.31 but since most of it is also relevant 
for more general modular abelian varieties considered in the present paper, we defer a discussion of 
the known results to the end of this introduction. 

We note that in the case of elliptic curves, the classical BSD conjecture (Conjecture [L2]) shares 
many of the same arithmetic quantities with the p-adic BSD conjecture (Conjecture II. 3^ : the main 
difference is that the regulator and L-series are replaced with p-adic analogues. In particular, the 
conjectures are equivalent if the rank is and p has good ordinary or nonsplit multiplicative reduction. 
Consequently, one might expect that a statement like Coni ecture 1 1 . 1 1 could be formulated and studied 
for a modular abelian variety A/Q associated to a newform /; this is the goal of the present paper. 

One of the main difficulties in formulating the conjecture is that the correct definition of the 
p-adic L-function of A involves some subtleties to get the precise normalization. By the general 
motivic framework outlined by Coates [10] , the p-adic L-series of A should interpolate the complex 
L-series L{A,s) at special values, so it seems plausible that it can be defined, similar to L(A, s), as 
the product of p-adic L-functions Lp{f'^,s) associated to the Galois conjugates of /. However, 
there is no obvious canonical choice for Lp{f'^, s), since picking a p-adic X-function associated to 
requires picking a Shimura period (see Theorem 12. 2p for each f^, that is, a complex number fi^^ 
such that L(/'^, is algebraic. In the case of elliptic curves, this is not an issue, since we can 

choose iljr to be the real period of the associated elliptic curve E. 

On the other hand, an extension of Coniecture ll.3l to modular abelian varieties should be equivalent 
to Conj ecture 1 1 . 1 1 in rank 0. Since the latter involves the real period associated to A, we normalize 
our Shimura period and hence our p-adic L- function, by requiring that the product of the 
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Shimura periods equals that we can do this is a consequence of a result we prove later (see 
Theorem [Q]). 

In this way we are able to explicitly construct, essentially generalizing the treatment in a 
p-adic i-function associated to A with the expected interpolation property, in the case where p 
is a prime of good ordinary reduction (sec (|2.10p and p.lip 'l. We keep the notation introduced 
prior to Conjecture 11.31 and extend it to the case of modular abclian varieties. We also define the 
p-adic multiplier ep{A) as follows: fix a prime p | p of the number field Kf generated by the Hecke 
eigenvalues of / and let a"' denote the unit root of — a{ap)x + p E {Kf)p[x], where a : Kf ^ C 
is an embedding. For a Galois conjugate of / define ep(/'^) = (1 — l/{a°'))'^ and define ep{A) to 
be the product of the p-adic multipliers Cpif) over all distinct Galois conjugates of /. 

We make the following p-adic BSD conjecture: 

Conjecture 1.4. Let A/Q he a modular abelian variety associated to a newform f and let p be a 
prime number such that A has good ordinary reduction at p. Then the Mordell- Weil rank r of A 
equals OTdT{Cp{A, T)) and 

\m{A/Q)\ ■ Reg^(^) • n„c„ 



where C*{A,0) is the leading coefficient of the p-adic L-series Cp(A,T). 

Note that the conjecture of Mazur-Tate-Tcitelbaum in the good ordinary case is a special case 
of Conjecture 11.41 Moreover, if the rank of A/Q is zero, then our conjecture is equivalent to the 
classical BSD conjecture due to the interpolation property 

Cp{A,0) = Lp{A,l) = ep{A) 



Most progress toward proving the Mazur-Tate-Teitelbaum conjecture uses Iwasawa theory and many 
results remain valid in our setup. Unfortunately, Iwasawa-theoretic results typically only assert 
equality up to a p-adic unit, whereas our Conjecture II .41 asserts full equality. We restrict to the good 
ordinary case from now on, although most results have supcrsingular or multiplicative analogues. 
See [43l §6,7] for a summary of such analogues in the elliptic curves case. 

Let X{A/oaQ) denote the Pontryagin dual of the p-Selmer group of A and let A be the completed 
group algebra Zp[[r]]. Using p-adic Hodge theory, Kato [H] has shown that X{A/ooQ) is a torsion 
A-module. Hence we can associate a characteristic series fA{T) e Zp[[T]], well-defined up to a factor 
in Zp[[T]]^ to X(A/ooQ). 

Conjecture 1.5. (Main conjecture of Iwasawa theory for abelian varieties with good ordinary re- 
duction) There exists an element u{T) € such that 

Cp{A,T) = fA{T)-u{T). 

If A = _E is an elliptic curve, then the main conjecture is known to be a theorem in many cases. If 
E has complex multiplication, then a proof is due to Rubin |35| . Many other cases have been proven, 
culminating in the work j40j of Skinner and Urban. See |43| §7] for an overview. The following 
unconditional result is due to Kato [18] : 

Theorem 1.6. (Kato) Let A he an elliptic curve. There is an integer m > such that fA{T) divides 
p"'Cp{A,T). 

The following result of Pcrrin-Riou [31] and Schneider [37] relates fA{T) to the right hand side 
of Ol: 
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Theorem 1.7. (Perrin-Riou, Schneider) The order of vanishing OTdT=o fA{T) is greater than or 
equal to the rank r of A/Q. Equality holds if and only if the p-adic height pairing on A is nondegen- 
erate and the p-primary part of the Shafarevich-Tate group of A is finite, in which case 

the leading coefficient of fA{T) has the same valuation as 

... |m(A/Q)(p)|.Reg^(A).n.c. 

\A{Q)tors\-\A^{Q)tors\ ' 

See for instance [37l Theorem 2'], noting that we have ordp{ep{A)) = 2ordp(7Vp), where A^j, is the 
number of Fp-rational points on the reduction of A over Vp. 
As a corollary of Theorem 11.61 and Theorem ll.7( we have that 

(1.3) ordT=o Cp{A, T) > otAt=o fA{T) >r = rank(A(Q)), 

if A is an elliptic curve, so one direction of the first part of Conjecture 11.31 (with p.3p suitably 
modified in the case of multiplicative reduction) is already known. Moreover, the second part of 
our Conjecture 11.41 is consistent with the Main Conjecture, since the latter implies that the leading 
coefficients of the p-adic L-series and the characteristic series have the same valuation. If Cp{E, 0) ^ 
0, then Conjecture 11.31 part (ii) is also known up to a rational factor (see [TB], [3T], and also the 
exposition in [431 §8]). The primes appearing in this factor can be determined explicitly using |40] . 
In the good ordinary case a similar result also holds when the p-adic analytic rank is 1, under an 
additional hypothesis; see [43l §9]. This follows from work of Perrin-Riou [32] and Kato |18| . 

Historically speaking, numerical evidence played a crucial role in the formulation of Conjecture [TT3] 
[24l §11.12]. Gathering evidence for Conjecture 11.41 would require two computations independent of 
the usual Birch and Swinnerton-Dyer conjecture: the computation of p-adic regulators of A, as well 
as the computation of special values of the p-adic L-function attached to A. We give algorithms to 
compute these quantities and provide the first numerical verification for Coniecture ll.ll bv considering 
the modular abelian varieties of dimension 2 and rank 2 in [13] and the Jacobian of a twist of Xo(31) 
of rank 4. 

Our aim is to give a self-contained discussion of the p-adic Birch and Swinnerton-Dyer conjecture 
for modular abelian varieties. To that end, we discuss both the theoretical and the algorithmic 
aspects of p-adic special values and p-adic regulators. This paper is structured as follows: in §2, 
we give a construction of the p-adic L-series attached to modular abelian varieties, making explicit 
certain aspects of [24|. This allows us to compute p-adic special values. In §3, we take a look at 
the p-adic regulator attached to an abelian variety, focusing on the case when the abelian variety is 
the Jacobian of a hyperelliptic curve. We begin by reviewing the work of Coleman and Gross [T^] , 
which gives the p-adic height pairing on Jacobians of curves in terms of local height pairings. We 
discuss the two types of local height pairings which arise and give an algorithm to compute p-adic 
heights, which allows us to compute p-adic regulators. In §4 and §5 we present the evidence for the 
conjecture in dimension 2. 
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2. p-ADIC L-FUNCTIONS ATTACHED TO MODULAR ABELIAN VARIETIES 

In this section, we construct the p-adic L-function attached to a ncwform / g S2(J'i(N)) with 
trivial Dirichlet character, making explicit a few aspects of [24l §1.10]. This depends on the choice 
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of a Shimura period as in Theorem 12.21 In order to pin down the p-adic L-function we want, we 
relate the Shimura periods of / and its Galois conjugates to the real period of the abelian variety A / 
attached to /. This leads to a definition of a p-adic L- function for Af which satisfies the expected 
interpolation property (j2.13p . Finally, we discuss how this p-adic L-function can be computed in 
practice. 

2.1. Periods. Let be a positive integer and let Xi{N) be the modular curve of level A^. The 
Jacobian Ji{N) of Xi{N) is an abelian variety over Q of dimension equal to the genus of Xi{N), 
which is equipped with an action of the Hccke algebra T. The space S := S2{Ti{N)) of cusp forms of 
weight 2 on ri(A^) with trivial Dirichlet character is a module over T. Let f{z) ~ X^^i ^ns'^™''' G S 
be a newform, let Kf be the totally real number field Q(. . . , a™, ■ . ■), and let // denote the annihilator 
AnnT(/) of / in T. Following Shimura [38], we have that the quotient 

Af = MN)/IfMN) 

is an abelian variety over Q of dimension g ~ [Kf : Q] which is equipped with a faithful action of 
T///. Moreover, A/ is an optimal quotient of Ji{N) in the sense that the kernel of Ji{N) Af is 
connected. For ease of notation, we will drop the subscript / and write A = Af. 

Remark 2.1. We assume that / has trivial character for convenience, because we need that Kf is 
totally real in order for Theorem 12.21 and Theorem 12.31 to hold precisely as stated. However, in the 
case where Kf is a CM-field, Shimura [33] has proved that a slightly modified version of Theorem 1 2. 2 1 
continues to hold. Using this, one can prove a result that is analogous to Theorem 12.31 for arbitrary 
newforms / S 52(ri(7V)). 

There is a complex- valued pairing ( , ) on S* x Hi{Xi{N),Z), given by integration: 

(ft., 7) = 27ri / h{z)dz. 

This pairing induces a natural T-module homomorphism 

$ : i7i(Xi(iV),Z) ^ HomclS", C), 

called the period mapping. 

Let G f he the set of embeddings a : Kf — > C . If ct G G/, we let J"' denote the conjugate 
of / by a. We denote the complex vector space generated by the Galois conjugates of / by S"/. 
Let <&/ : Hi{Xi{N),Z) — >■ Homc(5'/,C) be given by 4> composed with restriction to Sf. Then 
^ f{Hi{Xi{N),Z)) is a lattice in IIomc(5'/,C) and we have an isomorphism 

A{C)^Romc{SfX)/^f{Hi(X,{N),Z)). 

A choice of basis B of Sf induces an isomorphism Homc(S'/, C) = and B maps via $/ to a lattice 
Ai3 C Cf such that 

^(C) = CVAfj. 

For a basis B of Sf we let Ag (resp. A^) be the fixed points of A^ under complex conjugation (resp. 
under minus complex conjugation). 

We define the real period f2j^ (resp. the minus period f2^) of A as follows: Let uia be the puUback 
of a generator of the nontrivial global relative differential g-forms on the Neron model ^ of ^ over 
Spec(Z) to A. We call lja a Neron differential on A. Then we define 

/ M 

Ja{C)± 

where ^(C)^ denotes the set of points of ^(C) on which complex conjugation acts as multiplication 
by ±L 
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Let ;B be a Z-basis of the finitely generated free Z-module consisting of elements of Sf with Fourier 
coefficients in Z. Then we have (see [U §3.2]) 

= p-ca-vo1(a|), 

where ca is the Manin constant of A, defined, for instance, in [U §3.1] and p G C is if ± = — and I 
otherwise. It is known that ca is an integer and conjectured that it is always 1 (cf. [U §3.3]). There 
is no known algorithm to compute the Manin constant in general, which complicates much of what 
wc do below. The evidence that c^i = 1 is compelling, and we make the following: 

Running Hypothesis: We assume for the rest of this paper that ca = ^■ 

Let Hi(Xi(Af), Z)* denote the part of Hi{Xi{N),'L) fixed by complex conjugation (rcsp. minus 
complex conjugation). If w, z e C, then we write w ^ z ii w and z differ by a rational factor. We 
have that ^ 

- p • vol {Kb ) , 

where Ag is the lattice <^ f{Hi{Xi{N), ^)*)i ^ f is induced by the choice of basis B as above, and 
the rational factor is the number of components of j4(M) (rcsp. A{C)^). 
We denote the complex L-series of / by 

n>l 

If ■)/) is a Dirichlet character, we denote its Gauss sum by T{tp) and its conjugate character by -ip. We 
also let /0 denote the newform / twisted by ip and K^f, the field generated over Q by the values of ip- 

Theorem 2.2. (Shimura, Theorem 1]) For all a £ Gf there exist iljr^ e M. and QJ„ G i-M such 
that the following properties are satisfied: 
(i) We have 



r{z)dz± J r(z)dzj eKf 



for all r e Q. 
(ii) If ip is a Dirichlet character, then 

L{f^A) 



-r £ Kf ■ Kj,. 

(Hi) If Tp is a Dirichlet character, then 



We call a set {^l^^lo-eG/ as in Theorem l2.2l a set of Shimura periods for /. Note that the conditions 

of Theorem l2.2l do not determine the sets {^^^}aeGf Indeed, if {^^^^}ai£Gf satisfy the assertions of 

the theorem, then this also holds for {cr{b) ■ fl^„}a-(^Gf, where 6 € 

± 

A- 

Theorem 2.3. Let {flf^j^fzCj be any choice of Shimura periods as in Theorem \2.2[ Then we have 



According to Shimura [39l §2], the periods ftf,, are related to a certain period lattice, which gives 
us a way to compare them to the periods £7^ 



- n 
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Proof. Fix a Z-basis B = {hi, . . . , hg) of the free Z-module consisting of elements of Sf with Fourier 
coefficients in Z. Then there are bi, . . . ,bg £ K such that 



If (7 e G/, then 



and hence wc have 

(2.1) (.r,7) = E'^(^^)-(^-^) 

1=1 

for each 7 e i/i(Xi(7V), Z). 

Now fix some ordering ai,...,ag of G/ and let B = (fey) be the g x 17- matrix with entries 
bij = (Tj{bi) We have that B' = {f'^)aGGf forms a basis of the complex vector space Sf. Wc will 
compute vol (Ag,) in two different ways and the desired equality up to a rational number will fall 
out. 

First we express vol (A^,) in terms of fl^. Note that (|2.ip implies 

(2.2) p ■ vol (A|,) - p • vol (a^^) - p • I det(B)| • vol (Ab"^) | dct(S)| • rj^. 

By [Sni §2], we can also relate the lattice Ag' to the Shimura periods ^f^. Namely, we let Kf 
act on as follows: if a e Kf, then the action of a on is represented by the diagonal matrix 
diag((cri(a))i). We can decompose 

Ab'^Q^ {Kb- ® Q)+ © (Ab' ® Q)-, 

where (Ag/ ® Q)^ is the set of elements of Ag/ fixed by ± complex conjugation. Then we have 

(As'®(Q)± = A7r!±, 

where fi^ S is the vector whose z-th entry is fij,,; . Hence (Ag/ (X) Q)^ is a one-dimensional 
i^/- vector space. 

Since det(i?) 7^ 0, the elements bi, . . . ,bg of Kf form a basis for Kf over Q and therefore a basis 
of (Ag/ (g) Q)* as a Q- vector space is given by 



(a-i{bi)flp,, . . .,ag{b,)VL\ 
Hence we have 

(2.3) p-vol(A±) ~ |det(i?)|-n^^/.> 



since 



{Kb'®W = (A^ 

The proof of the theorem now follows from (|2.2p and (|2.3p . □ 

Remark 2.4. If L{f, 1) 7^ 0, then we can also argue as follows: We have 

L{A,\) 

— \ e Q 
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by Theorem 12.21 But on the other hand, the quotient is a rational number as well by [H 

Theorem 4.5]. Hence we get HcreG/ ^ 

Remark 2.5. In [46], Vatsal defines canonical Shimura periods associated to cuspforms. It would 
be interesting to determine whether his periods satisfy Theorem 12.31 

From now on, we fix some choice {^^„}aeG f such that 
(2.4) J] n%=n% 

If {^'jcrjcrGG/ is another set of Shimura periods satisfying (|2.4I) . then there is a unit b G 0Kf such 
that ^'^<, = cr(6)rij„ for all cr G G/. For our intended applications, this ambiguity is not serious, see 
Remark [in 

In order to compute {r2^„}o-eG/ we can find a Dirichlet character ■0 of sign ± such that L{f^, 1) 7^ 
and use equation (11) of [JU]- Alternatively, we can fix some nonzero element 7 G Hi{Xi{N), Z)[If]^ 
and define 

f}± =a(6).(r,7), 

for each a G Gf, where b G Kf is chosen to make (|2.4p hold. See 01] Ch. 10] for a description of 
how to compute the integration pairing in practice. 

As an application of Theorem 12.31 we can prove a relation between the real and minus period of 
A and the corresponding periods of A twisted by a Dirichlet character ip. 

Corollary 2.6. Let be a Dirichlet character such that L{f^, 1) 7^ 0. Then there exists rj^ G 
such that 

In particular, if tJj takes values in Q, then there exists rj^ G Q* such that 

Proof. This follows from Theorem [131 and Theorem [221 □ 

Remark 2.7. If ^ is the Jacobian of a hyperelliptic curve of genus at most 2 and -i/) is a quadratic 
Dirichlet character such that ip{N) 7^ 0, then one can show that the statement of Corollary 12.61 
holds without the assumption L(/^, 1) 7^ using quite concrete arguments. For elliptic curves and 
quadratic -0, CoroUarv 12.61 was already used in (Ml §11.11]. Note, however, that their claim that 
77^ G {1,2} is incorrect; see |30j . where the correct value of rj^ is determined in all cases. 

2.2. Modular symbols, measures, and the p-adic L-function of a newform. In this subsec- 
tion we define the p-adic L-function associated to /, following [H]. See also the treatment in 
The definitions for f"', where a Gf, are entirely analogous. 

Recall that we fixed a choice of Shimura periods fij^ above. The plus modular symbol map 
associated to / is the map 

' ^ = fi^)dz + f{z)dz^ 
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and the minus modular symbol map associated to / is the map 



J/ 



r H- > [?'] 



/ 



/ 



f{z)dz ■ 



f{z)dz 



Note that we have [0]+ ^ 



More generally, if m is a positive integer and is a, Dirichlet 
character modulo m, then j24[ §1.8] implies 



(2.5) 



E 

u mod m 



V(a) 



e K 



/■ 



Let p be a prime of good ordinary reduction for A. We fix, once and for all, a prime p of Kf 
lying above p. The modular symbol maps allow us to define two measures on which depend on 
the unit root of the polynomial h{x) :~ — QpX + p e {Kf)p[x], where {Kf)p is the completion of 
Kf at p. The construction of the p-adic L-function depends, in turn, on these measures. Since A is 
ordinary at p, the polynomial h has a unique unit root a G (if/)p, i.e., a root with ordp(Q;) = 0. 

Using the modular symbol maps [ ]j , we define two measures ^ on Z^ by 



M^,„(a + p"Zp) 



-i,ri+l 



P" 



For a continuous character x on ^p with values in Cp, we may integrate x against /^/,q,. Following 
[24l §1.13], we write x g Z^ as a;(a;) • (x) where w(x) is a (p — l)-st root of unity and {x) belongs to 
1 +pZp. The element w is known as the Teichmiiller character. 

We define the analytic p-adic L-function associated to / by 



Lp{f,s)= / (.t)-' ^d^ii{x) foraUseZp, 



(2.6) 

where by (x)'^^^ we mean expp((s — 1) • logp(a;) and exp^ and log^ are the p-adic exponential and 
logarithm, respectively. The function Lp(f, s) extends to a locally analytic function in s on the disc 
defined by |s — 1| < 1, as in the first proposition of [MJ §1.13]. 

Let ooG be the Galois group Gal(Q(/ipoo)/Q). The cyclotomic character k : — )• Z^ induces an 
isomorphism from ooG to Z^ that sends a topological generator 7 in ooG^^~^'' to a generator ^(7) of 
1 
T 

(2.7) Cp{f,T)= I (l + T)^^^^ dfi+Jx). 

Now for each n > 1, let P„(/, T) be the following polynomial: 



pZp . This identification allows us to give a series expansion of the p-adic L-function in terms of 
: ^(7)''^"'^ — 1. That is, we have 



p-i /p"-i-i 

(2.8) P„(/,T) = ^ Yl M;,„(^(«)(l+p)^'+p"Zp).(l + r)^- 

a=l y j=0 

We have that (|2.8p gives us a Riemann sum for the integral (|2.7p , by summing over residue classes 
mod p"; in other words: 

Proposition 2.8. We have that the p-adic limit of these polynomials is the p-adic L-series: 

hm PM,T)^Cpif,T). 
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(2.9) Cpif^O) - Lp(/, 1) = / dM7,„ = ^pif) ■ mj = epif) 



This convergence is coefficient-by-coefHcient, in the sense that if P„(/, T) = J2j^n.jT^ and 

hm a„j = aj. 

Proof. This is a straightforward generahzation of |431 Proposition 3.1]. The upper bounds we obtain 
are the same as the upper bounds in the proof of |431 Proposition 3.1], which enables us to compute 
the p-adic L-serics to any desired precision. □ 

We define the p-adic multiplier ep{f) by 

The p-adic L-series of / satisfies an interpolation property with respect to the complex L-series 
of / [Ml §1-14]: 

-v J 

2.3. p-adic L-function associated to A. The p-adic L-scrics Cp{f,s) associated to / that we 
constructed in the previous subsection depends on the Shimura periods f2j^„ and on the prime p. 
In the present section we define a p-adic i-function associated to the abelian variety A which is 
independent of the choices of the Shimura periods (provided they satisfy (|2.4p ) and of p. 
The abelian variety A has an associated complex L-serics, given by 

L{A,s)^ n L{r,s), 

which can be extended analytically to the whole complex plane. 
We define the p-adic L-function associated to A by 

(2.10) ip(A,s)= H Lpir,s) 

aeGf 

for s G Zp. 

Remark 2.9. Since we require our Shimura periods {Cl^^}„^Cf to satisfy (|2.4p . the p-adic L- 
function Lp(A, s) does not depend on the choice of period for each f , although the individual p-adic 
L-functions Lp{f"',s) do. 

Furthermore, we define 

(2.11) CpiA,T)= H Cpir.T) 



creGf 



and the p-adic multiplier of A by 



For ?■ e Q we set 

(2.12) [r]± := [] 



aeGf 



Corollary 2.10. Let {fl^^}a-eGf be a set of Shimura periods satisfying (|2.4|1 . Then (|2.9|) immediately 
implies 

(2.13) Cp{A, 0) = Lp{A, 1) = ep{A) ■ [0]+ = ep{A) ■ 
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2.4. Quadratic twists and normalization. Modular symbols can be computed up to a rational 
multiple purely algebraically (cf. |41| ) using (mostly sparse) linear algebra over fields. Computing 
the exact modular symbol (not just up to a rational factor) requires doing linear algebra over Z. 
which is much slower. In this section we describe a method to determine the correct normalization 
of the modular symbol map by using special values of quadratic twists, which is potentially much 
faster than using linear algebra over Z. 

In order to find the correct normalization, we use the fact that the p-adic i-series associated to 
A interpolates the Hasse-Weil L-function L{A, s) associated to A at special values. Algorithms for 
the computation of L{A, 1) and 0+ are discussed in [3 [T31 SI] . So if L{A, 1) ^ 0, then we 

A 

can find the correct normalization factor S'^ for £.p{A,T) by computing [0]^ and comparing it to 
■ Note that the quotient '^[^''"^ can be computed purely algebraically as a certain lattice index, 

without computing either of the real numbers L{A, 1) or 

In order to discuss the strategy for the case L{A, 1) = 0, we begin by considering modular symbols 
associated to quadratic twists of /. Let _D be a fundamental discriminant of a quadratic number 
field such that gcd{pN, D) = 1 and let ^ denote the Dirichlet character associated to Q{\/D). We 
assume that L{f^, 1) 7^ 0. Using Corollarv l2.6l we see that there exists % € K'j: such that 

Qsign(D) _ 0+ 

A computation analogous to (331 §3.7] yields 

o-eG/ crGG/ «=1 

Therefore we can compute the product of the plus modular symbols for and its conjugates in 
terms of modular symbols for / and its conjugates. The same holds for the p-adic L-function of the 
twist A^ of A by '0- 

Now suppose that Cp{A, 0) — and we want to find the correct normalization factor for [ ]j . We 
can use that for a fundamental discriminant D with Dirichlet character ip, the modular symbols [ ]y 
and [ J/^"'"^'' are related by (|2.14p . Hence the same normalization factor (5+ will yield the correct 
value [r]^ for all > such that gcd{pN, D) = 1. 

We can compute by finding a fundamental discriminant D > such that gcd{pN, D) = 1 and 
such that A^ has analytic rank over Q, where tp is the quadratic character associated to Q(\/15), 
and comparing [0]^^ to 

r]^ ■ L{A^, 1) 

An analogous approach can be used to find the correct normalization factor for the minus modular 
symbol. It follows from [S] that in both cases a fundamental discriminant D as above always exists. 

Remark 2.11. Suppose that A is the Jacobian of a hyperelliptic curve X/Q of genus g given by an 
equation 

+ hi{x)y = fi{x) 

which is minimal in the sense of [2^. Often it is not necessary to compute il^®" (or even n-X) to 
compute r]^. Let + h2{x)y = f2i^) be a minimal equation for X^,, and consider the differentials 
uji ~ 2y+hi{x) '^^ ^^'^ ^'i ~ 2y+h2{x) frequently happens that (wi, . . . , Wg) is a basis of 

the integral 1-forms on A and (wj^, . . . ^uj'g) is also a basis for the integral 1-forms on A^. In that 
case we always have 77^ e {±1}; this follows from [13l §3.5]. It is easy to determine the sign using a 



u 
D 
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straightforward generalization of [l9l §1.3]. More generally, a similar approach can also be used to 
compute rjip directly if we know how to express a basis for the integral 1-forms in terms of cji, . . . , Wg. 

2.5. The algorithm. We implemented the following algorithm for the computation of the p-adic 
L-series of A in Sage [42] . 

Algorithm 2.12 (p-adic L-series). 

Input: Good ordinary prime p, A modular abelian variety attached to newform /, precision n. 
Output: nth approximation to the p-adic L-series Cp{A,T). 

(1) Fix a prime p of the field Kf generated by the Hecke eigenvalues of / lying above p and 
compute the unit root a of h{x) G {Kf)p[x]. 

(2) Find a fundamental discriminant D > such that gcd{pN, D) ^ 1 and A^ has analytic 
rank over Q, where -0 is the quadratic character associated to Q{\/D). 

(3) Compute rj.^. 

(4) For each a G Gf, define measures ^j,, ^. 



Remark 2.13. Note that Step (7) of Algorithm l2.12] is exponential inp; sec the following subsection 
for an alternative method. 

2.6. Overconvergent modular symbols. Here we outline an alternative method for Step (7) of 
Algorithm l2.12l This method is due to Pollack and Stevens [34] and has running time polynomial in 
p and in the desired number of digits of precision. 

The idea is to use Stevens's overconvergent modular symbols; these are constructed using certain p- 
adic distributions, and they can be specialized to classical modular symbols. More precisely, any clas- 
sical modular cigensymbol can be lifted uniquely to an overconvergent modular Hecke-eigcnsymbol, 
which can be approximated using finite data. Note that in order to do this, we first have to p-stabilize 
the symbol to a symbol for ro{Np) which is an eigensymbol away from p. 

The plus modular symbol we start with is only determined up to multiplication by a scalar, so 
the corresponding overconvergent eigenlift is also only determined up to multiplication by a scalar. 
Hence we cannot dispense with Steps (2), (3), (5) and (6) of Algorithm 12. 121 

Once this desired lift has been computed, writing down the p-adic L-series associated to the 
modular symbol and its quadratic twists by ip for suitable D is rather easy, cf. |34[ §9]. Together 
with David Roe and Robert Pollack, we have implemented Algorithm 12.121 with Step (7) replaced by 
the algorithm from [34] in Sage as well, building on an implementation due to Pollack. 

3. The p-ADIC height pairing of COLEMAN-GrOSS and p-ADIC REGULATORS 

We now shift our attention to the remaining p-adic quantity appearing in Conjecture II. 4i the p- 
adic regulator. To discuss p-adic regulators, we begin, in this section, by describing one construction 
of the global p-adic height pairing relevant to our setting. We give an algorithm to compute the 
height pairing in the case when the abelian variety A is the Jacobian of a hyperelliptic curve and 
show how we use it to compute p-adic regulators. 

Let A be an abelian variety defined over a number field K and let denote the dual abelian 
variety to A. There are several definitions of p-adic height pairings on abelian varieties in the 
literature. Schneider [36j and Mazur-Tate |23j were the first to construct the p-adic height pairing on 
abelian varieties defined over number fields. This was extended to motives by Nekovaf [29]. There 





(6) 

(7) 

(8) 




using (I2T3]) . 
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are also more specialized definitions: in the case when dimyl = 1, K ~ Q, and p is a prime of 
good, ordinary reduction, Mazur, Stein, and Tate |22| gave an explicit formula for the p-adic height 
which relics on an understanding of the p-adic sigma function. When A is the Jacobian of a curve, 
Coleman and Gross [12] described the p-adic height pairing on ^ as a sum of local height pairings. 
Note that in the range where all of these constructions apply, they are known to be equivalent by 
the work of Coleman [llj and Besser [6] (where the equivalence is possibly up to sign, e.g., in the 
supersingular case). For all of these definitions, the p-adic height pairing is known to be bilinear 
and, in the principally polarized case, symmetric. 

Let p be a prime number such that A has good ordinary reduction at all primes of K above p. 
We denote the p-adic height pairing by 

h: Ax A"" — >Qp 

{P,Q)^h[P,Q). 

Now we want to define the p-adic regulator with respect to h. In the literature, one usually defines 
this quantity as the determinant of the height pairing matrix with respect to set of generators of 
the free part of A{K) and A"^ {K), respectively. This, however, is only well-defined up to sign. Since 
results in Iwasawa theory are typically only up to a p-adic unit, this is usually not a serious problem, 
but in order to state Conjecture II. 4[ we need a canonical well-defined p-adic regulator. First note 
that \i (f) : A ^ is an isogeny and Pi, . . . , generate the free part of A{K), then 



Keg^{A/K) :=det(^(/i(P„0(P,)))^^^. 
does not depend on the choice of generators Pi, . . . , P^. 

Lemma 3.1. Let c,c' € Pic'^(A) and let 4'c,(t>c' '■ A — > A"^ be the corresponding isogenies. If 
Pi, . . . , Pr is a set of generators for the free part of A(K), then we have 

^ . Reg^^ (A/K) = ] ,,,,,,,, Reg. , {A/K) . 



[A^{K):MAiK))] ' ' [A^{K):cj^AA{K))] 

Proof We set m = [A^ {K) : (j)c{A{K))] and m' = [A^ {K) : (jy^'iAiK))]. Let Qi, . . . , denote a set 
of generators of the free part of A'^{K) such that 

det ((/i(P„Q,)),,,.) ^^Reg^M/K). 
There are integers mij, m,^^ e Z such that 



= (Pi) = '^rriijQj and (/)c'(Pi) = 

j=i i=i 



''3- 



Since dct (j^rriij) = m and det ^(/i(Pi, Qj))- ^-^ = ±:^Hcg^^{A/ K), it suffices to show that 

det ^(m-j). = m' (and not —m'). 

To show this, we use classical Neron-Tate heights. The Neron-Tate regulator is defined as 

Rcg(A//0 = |det(((P„Q,)p)^^^.)|, 

where the bilinear pairing {Pi,Qj)-p = h-p{Pi,Qj) is defined in terms of the Neron-Tate height h-p 
with respect to the Poincare bundle on A x A^. See [111 Remark F.4.1]. On the other hand, it is 
easy to see that 

{p,,MPj))v^{n,Pj)c, 
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where ( , )c is the Neron-Tate height pairing on A with respect to c, and similarly for c'. This implies 

that 

< det [m,P,)c),^,) = det(Af) • det (((P., Q,)p).,^-) . 
Since dct(M) = m > and Rcg(A/A') — Rcg^(A//C), we have 



Reg(A/A')=det(^((P„Q,)p),^^.^ 
and hence 

Keg(A/AO = ^dct (((P.0.(P,)).).,) = ^dct . 

Therefore we find det(M') = m! . □ 

Definition 3.2. Let A be an abelian variety defined over a number field K and let denote 
its dual. Let Pi, . . . ,Pr be a set of generators of A{K)/A{K) tors- Fix some c £ Pic°(A) and let 
(j)c '■ A ^ denote the corresponding isogcny. The p-adic regulator of A, denoted Rcgp(yl), is 
defined by 

^^s^(^) ^- [Av(,,^)/^^(^(^-))] Reg,.(A/j^) 

The p-adic regulator is well-defined by Lemma [01 It has been conjectured by Schneider [35] that 
the p-adic height pairing is nondegenerate, however, in contrast to the classical case of Neron-Tate 
heights, this is not known in general. 

Among the aforementioned definitions of the p-adic height pairing, the Coleman-Gross construc- 
tion of the p-adic height pairing is fairly explicit in nature, and for that reason, lends itself nicely to 
computation. Thus we take it as our working definition of the p-adic height. We start by giving a 
brief overview of the work of Coleman and Gross. 

Suppose X/K is a curve defined over a number field K, with good reduction at primes above p. 
To define the p-adic height pairing 

h : Div°(X) X Div°(X) ^ Qp , 

where Div°(X) denotes the divisors on X of degree zero, one needs the following data: 

• A "global log"- a continuous idele class character £ : A*j^/K* — >■ Qp . 

• For each v \ p a. choice of a subspace C iJjfj((X (g) K^)/K^) complementary to the space 
of holomorphic forms. 

We require that the local characters £y induced by £, for v \ p, are ramified in the sense that they do 
not vanish on the units in Ky. From £ one deduces the following data: 

• For any place v \ p we have £v{&k ) = for continuity reasons, which implies that £y is 
completely determined by the number ^^(7r.„), where 7r„ is any uniformizer in Ky. 

• For any place v \ p we can decompose £y as a composition 

(3.1) 0t ^Qp 




Ky 

where ty is a Qp-linear map. Since we assume that £y is ramified it is then possible to extend 
log„ to log„ : K* — >• Ky in such a way that the diagram remains commutative. 
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We will later need to choose a branch of the p-adic logarithm, since the Coleman integral of a 
form with residue depends on such a choice. We will fix this choice for the computation of the local 
height pairing to be the one determined above. 

Let us now describe the p-adic height pairing h{D,E) for a pair of degree zero divisors D and E 
with disjoint support. The height pairing is a sum of local terms 

hiD,E) = Y,hv{D,E) 

V 

over all finite places v. The local terms depend only on the completion at v of K. Thus, let Ky be 
the completion of if at a place v, with valuation ring O, uniformizer tt and let ky = O/irO be the 
residue field, with order q. Let C denote the curve X over the local field Ky- We shall assume that 
C has a /Cu-rational point and that C has good reduction at tt. 

Let X '■ K* — > Qp be the local component of £, which is a continuous homomorphism. 

Proposition 3.3. //char ky ^ p, there exists a unique function {D, E) defined for all D, E £ Div'^(C) 
of disjoint support that is continuous, symmetric, hi- additive, takes values in Qp, and satisfies 

(3.2) {if),E)=xif{E)) 
forfeKy{C)*. 

Proof See [H Prop 1.2]. □ 

We will discuss how to compute this in practice in Section 13.11 

3.1. Computing p-adic heights away from p. Wc keep the notation of the previous section, 
but assume, in addition, that X is hyperelliptic of genus g, given by an equation = fi^), where 
f £ 0[x] is separable. Let u be a fixed non-archimedean place of K not dividing p. 

The arithmetic geometry needed in the present section can be found in [5T1 Chapters 8,9]. We 
fix a proper regular model C oi C = X Xk Ky over Spec(C') with special fiber Cy. If £) is a prime 
divisor on C, then we let D denote the Zariski closure of £> on C and wc extend this to all of Div(C) 
by linearity. It was shown by Hriljac [17] that if £> e Div*'(C), then there exists a vertical Q-divisor 
^{D) on C such that the intersection multiplicity oi D + ^{D) with any irreducible components of 
Cy is trivial. 

li D, E Div°(C) have disjoint support, then according to [12l Prop 1.2] the local height pairing 
between D and E at v is given by 

(3.3) hy{D,E)=ey{TTy)-iy (D + ^ (D) ,E) , 

where iy denotes the (rational- valued) intersection pairing on C. This does not depend on the choice 
of $(£») or of C. 

As in [26], the following steps are sufficient to compute the local p-adic height pairing at v. 

(1) Compute a desingularization C in the strong sense of the Zariski closure C of C over Spec(C'); 

(2) Compute iy (D^E) ; 

(3) Compute iy {<^>{D),E). 

These steps are dealt with in detail and greater generality in [26]. For the convenience of the reader, 
we provide a brief summary in the present case of hyperelliptic curves. 

Step (1) can be done using a desingularization algorithm implemented by Steve Donnelly in Magma 
[Hj. Recall that a desingularization C of the Zariski closure C of C over Spec(C') in the strong sense 
is a proper regular model of C over Spec(O) such that there exists a morphism ^ : C — > C that is 
an isomorphism above regular points of C. From now on we will assume that our model C is of this 
type, as this property is needed in order for some of the other steps to work. See [26l §4.3] for details. 
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For Step (2), we write our divisors D and E as differences of effective divisors 

D^Di-D2, E = Ei~E2. 

By bilincarity of the intersection pairing it suffices to discuss the computation of iu (-Di, i^i). 

For now we assume tliat tlie points on Cy where Di and Ei intersect aii fie on a singie affine 
piece of C. Suppose that = Spec(C'[a::i, . . . , Xn]/J) for some ideai J and that I-^ (resp. I-j^) 
represents Di (resp. Ei) on C. Then we have 

(3.4) (d;,E;) = iengthcj^g . . . , x„]/J + + " 

The computation of the right hand side of p.4[) can be reduced to (essentiaiiy) the computation of 
Grobner bases over Spec(C'), cf. [26l Aigorithm 1]. 

In order to find the representing ideals Ijy^ and /-g^ the strategy is to first find representing ideals 
for the Zariski closures of Di and Ei on C and lift these to C through the blow-up process. We can 
guarantee that the intersection of these closures has support only in one of the two standard affine 
pieces of C by decomposing Di and Ei into prime divisors over a finite extension M of K^; in the 
present case of hyperelliptic curves this is possible using factorisation of univariate polynomials over 
M as described in l26l §5.3]. 

So it remains to discuss how to represent Zariski closures of prime divisors on C on the affine 
piece 

r = Spec {0[x,y]/{y'-f{x))). 
If Di = X]iLi(-P'i) is reduced, then we can use a representing ideal 

{a{x),y - b{x)), 

where a{x) £ 0[x] has roots x{Pi), . . . ,x{Pd) and does not vanish modulo 7r„ and b{x) £ 0[x] does 
not vanish modulo 7r„ and satisfies y{Pi) = b(x{Pi)) for i e {1, . . . , d}. This is commonly referred to 
as Mumford representation, see [27l 3.19]. In particular, if Di = (-Pi), where Pi G C{Ky). then we 
can take the ideal 

{x~x{Pi),y-y{Pi)). 

The other case we have to consider is the case Di = (Pi) + (P^ ), where Pi is defined over an 
extension of Ky of degree at most 2 and P[^ is the image of Pi under the hyperelliptic involution. 
Then we can simply use the ideal 

{x^x{Pi)). 

For Step (3) we refer to [26l §4.5]. In brief, we first compute the intersection matrix M of Cy 
and its Moore-Penrose pseudoinverse M+. Suppose that the special fiber Cy is given by X^I^o ^i^i' 
where Fq, . . . are the irreducible components of Cy. We also need the vectors s{D) and s{E) of 
intersection multiplicities, where 

s{D) = (no • iy{D,ro), ■■.,nra- iy{D,T^)) ; 

and s{E) is defined similarly. These can be computed using the techniques introduced in Step (2) 
above. Then we have 

iymD),E) = -s{Ef ■ M+ ■ s{D). 

We have not discussed how we can compute a finite set U of places of K such that we have 
hy{D, E) ~ for all v ^ U. This is discussed in [26l §4.2, §5.2]. Here we only mention that it suffices 
to compute U containing all bad places (that is, all places v such that ord„(2 • disc(/)) > 0) and all 
places V such that D and E have nontrivial common support modulo tt^, . The latter can be computed 
as follows, where D ^ Di — D2 and E ^ Ei — E2 are as above. 
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We only discuss the computation of all v such that Di and Ei have nontrivial common support 
modulo TTi,. Let Idi and Ie^ denote representing ideals of the Zariski closures of Di and Ei on the 
affine piece 

of the Zariski closure of X over Spec(0/f ). We assume that Ok is Euclidean; the general case can 
be reduced to this situation using a straightforward trick discussed in [211 §4.2]. If S is a Grobner 
basis of the ideal 

[y^ - f{x)) + lD^+lE,. 

over Ok, then B contains a unique element q G Ok (cf. [Ml Lemma 4.3]). Factoring [q) yields 
a set of places containing all places such that Di and Ei have nontrivial common support on the 
reduction of the affine piece given by = f(x). Repeating this process for the other standard affine 
piece ~ x^^^"^ f{\/x) yields all places v such that Di and Ei intersect nontrivially modulo iTy. 

3.2. Computing p-adic heights above p. We now describe the local contribution at a place v \ p. 

Definition 3.4. Let D,E £ Div°(C) have disjoint support. The local height pairing at a place v \ p 
is given by the following Coleman integral: 



h,{D,E) :=t, (^J^ 



E 

where is the trace map determined by the decomposition of (see (j3.1[) ). and lod is a, differential 
associated to D. 

We start by reviewing the construction of lod- Let T{Kv) denote the subgroup of differentials 
on C of the third kind. We are interested in a particular subgroup of T{Ky) whose elements are 
the logarithmic differentials, i.e., those of the form ^ for / G K^{C)* . We denote this subgroup 

as Ti{Ky). Letting H^^{C / Ky) denote the space of holomorphic differentials and A = Pic"(C), we 
have the short exact sequence 

H\^{C/Ky) T{Ky)/Ti{Ky) A{Ky) 0. 

This sequence has a natural identification with the A'„-rational points of an exact sequence of 
commutative algebraic groups over Ky: 

H^^iC/Ky) ^ A^O, 

where U is the universal extension of A by a vector group and H^^{C / Ky) = G£. 

Now as Ky is p-adic, we will make use of the fact that we have a logarithmic homomorphism 
defined on an open subgroup of the points of any commutative p-adic Lie group, G, to the points of 
its Lie algebra Lie(G). When G = U oi A, the open subgroup on which the logarithm converges has 
finite index, so the homomorphism can be uniquely extended to the entire group. We denote this 
extension as log^ or log^, respectively. Since the logarithm is functorial and equal to the identity 
on H^^{C/Ky), we have the following: 

Proposition 3.5. There is a canonical homomorphism 

* : T{Ky)/Ti{Ky) Hl^{C/Ky) 

which is the identity on differentials of the first kind and makes the following diagram commute: 

^ h]:1{G/K,) ^ U{Ky) ^ A{Ky) ^ 

*=log„ log^ 

rl.O 



H'^^{C/Ky) ^ H\^[ClKy) ^ H'iC Oc/kJ 
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Note that the map ^' takes a differential of the third kind on C to a differential of the second kind 
modulo exact differentials, sending log differentials to 0. It can be extended to a linear map from 
the iir„-vector space of all differentials on C /Ky to H\^{C / Ky) by writing an arbitrary differential 
V as a. linear combination v = 'Yl, o^ifJ-i + 7j where fii is of the third kind, S Ky, and 7 is of the 
second kind on C. We then define ^{i^) "^^'(/ii) + [7]. 

Now recall that we have at our disposal the complementary subspace W = Wy. It allows us to 
isolate a canonical form ujjj with residue divisor D as follows: 

Definition 3.6. Let D e Div°(C). Then iojj is the unique form of the third kind satisfying 

rcs(a;£,) = D, e W. 

3.3. Computing the global p-adic height. Using the material in Sections 13.11 and 13.21 we may 

now give an algorithm to compute global p-adic heights on Jacobians of hyperclliptic curves. 

Algorithm 3.7 (Global p-adic height pairing). 
Input: 

• Genus g hyperclliptic curve X over Q of the form = f{x), with / separable, deg / = 2g + 1, 

• Prime p of good ordinary reduction for X, 

• Divisors D,E £ Div°(X) with disjoint support. 

Output: 

• Global p-adic height pairing h(D, E) = J^v ^v{D, E) 
Algorithm: 

(1) Heights away from p 

(a) Find relevant places. Compute the following set U of non-archimedcan places of K: 

U = {v : V \ 2 disc(/)} U {v : supp(£' mod v) n supp(i? mod v) ^ 0}. 

(b) Local computations. For each v €U go through the following steps. 

(i) Regular models. Compute a dcsingularization C of the Zariski closure of A 
Ky over O ^ Oy in the strong sense. 

(ii) Ideal representatives. Write D = Di — D2 and E = Ei — E2, where Di, Ej 
are effective and find representatives Ijj- and /-jr- of the Zariski closures Di, Ej of 
Di, Ej, respectively, on an affine piece of C containing supp(£') n supp(£'). 

(iii) Horizontal data. Compute iy{Di,Ej) for i,j G {1,2} using (|3.4|) . 

(iv) Vertical data. Compute the Moore-Penrose pseudoinverse of the intersec- 
tion matrix of the special fiber Cy = J^iLo niVi and for H € {D, E}, the vectors 

s{H) = {no ■ iy{H,ro), ...,n„i- iy{H,r,n)) ■ 

(v) Local height pairing. Let ky be the residue field at v and set 

hy{D,E)^ {^{E) ■ M+ ■ s{Df -Y,^v{Dl,'E^)^ • log(#A:,). 

(c) Global height pairing away from p. Compute J2y]p^v{D, E). 

(2) Height above p 

(a) From D to ljjj. Choose ui a differential of the third kind with Res(w) = D and compute 
log(u;) = \I'(w) for oj. Using the decomposition 

H'dR{C/Ky) ^ Hl'^{C/Ky) © W, 
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write 

log(w) = ?7 + \og{ujD), 

where 77 is holoniorplric, and log(cj£)) G W. Now, using this holomorphic component 77, 
we have 

ujo := Lo — rj. 

(b) Coleman integration 

(i) ...of a holomorphic differential. Compute 77, as in [4]. 

(ii) ...of a meromorphic differential. Let be a p-power lift of Frobenius and set 
a := (j)*uj — puj. Write E = Ei, where Ei = (Ri) — (Si) for points Ri, Si on X. 
Then for /3j a differential with residue divisor Ei, we compute 

. ^ _1_ U + 1: Rc, (a/ft) - £^ . - , 

as described in [3]. 

(c) Height pairing above p. Subtract the integrals to recover the pairing at p: 

hp{D,E)^ / a;^, = / cj - / 77. 

Je Je Je 

(3) Global p-adic height pairing. Return the sum of 1(c) and 2(c). 

Remark 3.8. Note that our current implementation of Algorithm 13.71 further assumes, in Step 
2(b)(ii), that R,,S, G C{Qp). 

3.4. Computing the p-adic regulator. In this section, we explain how we use Algorithm 13.71 in 
practice, to compute the p-adic regulator of A/Q, where ^ is a Jacobian surface of a curve X/Q of 
genus 2 and p is a prime such that A has good ordinary reduction at p. 

Suppose that P, Q G A{Q) are distinct and that we want to compute the p-adic height pairing 
of P and Q. If we can find representatives Di = (Pi) — (P2) and D2 = (Qi) — {Q2) of P and Q, 
respectively, where Pi,P2,Qi,Q2 G X{Q) are aU distinct, then we can simply apply Algorithm 13.71 
directly to compute hp(Di, D2). 

However, in many situations, it is not possible to find representatives of P and Q whose sup- 
port consists of Q-rational points. We can still compute the p-adic height pairing if we can find 
representatives Di,D2 G Div'^(Ar)(Q) of P and Q, respectively, such that 

Di XQ Qp = (Pi) + (P2) - (P) - (P-), 

D2 XQ Qp = (Qi) + {Q2) - (S) - (5-), 

where Pi, P2, Qi, Q2, P, G C{Qp) are pairwise distinct. It is explained in Section [QI how to find 
ideal representations of the positive and negative parts of Di and I?2, respectively, which is all we 
need to compute the local height pairings away from p. The latter is very similar to the computation 
of canonical real- valued height pairings on Jacobians of hyperelliptic curves as discussed in [2^ and 
has been implemented in Magma. 

If we want to compute the p-adic height pairing of P with itself, then we simply compute minus 
the p-adic height pairing of P with — P using the techniques discussed above. 
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Remark 3.9. In principle this algorithm can be generalized immediately to hyperelliptic curves 
over Q of arbitrary genus, with minor subtleties if the genus is odd and the degree is even. See the 
discussion in |26j . However, to find the p-adic regulator precisely (and not only up to a rational 
square), we first need a set of generators for A{Q)/A{Q)tors- Given generators of a finite index 
subgroup of A{Q)/A{Q)tors, a general algorithm for this computation is presented by StoU in [44], 
but currently it is only feasible if the rank is zero or X has genus 2 (using current work in progress due 
to StoU this can be extended to genus 3). Hence we can only hope to compute the p-adic regulator 
up to a rational square in general. 

Remark 3.10. It is useful to note that in order to compute the p-adic regulator in the genus 2 
case, we need not work with generators of the free part of the Mordell-Weil group itself. Indeed, 
if we have a set of points generating a subgroup G of finite index, then we can easily compute the 
index of G by computing and comparing the real- valued regulators of G and of ^(Q)/^(Q)tors, see 
Remark 13.91 Since the p-adic height is quadratic, it suffices to compute the p-adic regulator of G in 
order to deduce the p-adic regulator. This helps in finding points with representatives of the required 
form. 

Remark 3.11. Suppose that X is hyperelliptic and defined over a number field K. Then we can still 
compute the p-adic regulator exactly as above if all completions Ky at places v \ p satisfy Ki, ^ Qp. 

Finally, following j43j , we note that the p-adic regulator has a natural normalization from Iwasawa 
theory, coming from the choice of topological generator 7. This is done so that the global p-adic 
height depends only on the choice of isomorphism F — > Zp, instead of on the Zp-extension. This 
normalization is carried out by dividing h[Pi,Pj) by logp(K(7)), or, alternatively, since the p-adic 
regulator involves a basis of dimension r, by taking 



Regp(A) 



Reg^(^) = 



4. Evidence for rank 2 Jacobians of genus 2 curves 

As we now have algorithms to compute the p-adic regulator and p-adic L-series, we proceed to 
verify Conjecture 11.41 for specific abelian varieties, using BSD data from [T31 Table 2]. We take as 
our list of candidate modular abelian varieties A those appearing in [13] of rank 2. 

We were able to find generators represented by divisors whose support consists only of Q-rational 
points for each of the rank 2 Jacobians taken from [13] except for the one with level 167. Hence we 
used the easier first approach outlined in Section r3.4l to compute the p-adic regulator for all Jacobians 
except for the one associated to level 167. 

In order to compute the special values £p(^, 0), we used an implementation in Sage of the al- 
gorithm outlined in Section 12.61 Previously, we had used Algorithm 13.71 obtaining fewer digits of 
precision. The results agreed up to the precision obtained using the latter. 

Table [4T] is taken from [13] and contains minimal models (in the sense of [19]) for each curve from 
[IB] whose Jacobian variety has Mordell-Wcil rank 2 over Q, as well as the corresponding level A^. 

The data presented in this section proves the following: 

Theorem 4.1. Assume that Jor the Jacobians oj all curves in Table [^TT] the Shafarevich-Tate group 
over Q is 2-torsion. Then Conjecture \1.4\ is satisfied up to the respective precision specified in the 
tables below for the Jacobians of all curves in Table \4^.1\ at all good ordinary p < 100 satisfying the 
hypotheses of our algorithms. 

Remark 4.2. The assertion that all Shafarevich-Tate groups are 2-torsion for these abelian varieties 
follows from the classical conjecture of Birch and Swinnerton-Dyer by [T5] . 



p-ADIC BSD FOR MODULAR ABELIAN VARIETIES 



21 



For our algorithms, we take the integral models in |T3l Table 1] and do a change of coordinates 
to obtain the corresponding models of the form = f{x). We record both models in Table |4J] 



N 


\cl \ 1 


, h(x)) for integral model 




/(x) for 


y2 /(x) model 




67 


[x^ — 


X, X"^ + X + 1) 




x'' 


+ 4x^ 


+ 2x4 ^ _|_ „ 


2x + 1 


73 




- 2x^ + x,x^ + x'^ + 1) 




x^ 


- 2x5 


+ x4 - 6x3 ^ 2x2 ^ 


4x + 1 


85 




x^ + ix^ - 2x + 1, x^ + x^ - 


hx) 


x^ 


+ 2x5 


+ 7x4 _^ g^3 _|_ ;l3^2 


-8x + 4 


93 


(-2a; 


^ +x'^ + x^,x^ + x^ + 1) 




x*^ 


- 6x5 


+ 5x4 _^ g^3 _|„ 2x2 


hi 


103 


\x^ + 


x'^.x^ +x2 + 1) 




x*^ 


+ 6x5 


+ 5x4 _|_ + 2x2 


hi 


107 




X^ — X — 1, x^ + x^ + 1) 




x^ 


+ 2x5 


+ 5x4 ^ 2x3 „ 2a;2 _ 


-4x- 3 


115 


(2a;3 


+ X^ + X, X^ + X + 1) 




x^ 


+ 2x4 


+ 10x3 + 5x2 g^T 


hi 


125,A 


\x^ + 


2x"* + 2x^ + x^ - X - 1, x^ - 


Fx + 1) 


x^ 


+ 4x5 


+ 10x4 + 10x3 + 5x 


2 - 2x - 3 


133,B 


{—x^ 


+ x'^ - 2x3 2x2 _ 2x, x^ + 


X2 + 1) 


x^ 


- 2x5 


+ 5x4 _ g2,3 ^ ;l0x2 


- 8x + 1 


147 


{x^ + 


2x'' + x^ + x^ + 1, x^ + x^ - 


hx) 


x^ 


+ 6x5 


+ 11x4 + 6x3 + 5x2 


+ 4 


161 


{x^ + 


4x2 + 4x + l,x3 + X + 1) 




x^ 


+ 2x4 


+ 6x3 ^ 27x2 + 18x 


+ 5 


165 


\x^ + 


2x'' + 3x3 + ^2 - 3x, x^ + X 


2 + x) 


x^ 


+ 5x4 


- 168x3 + 1584x2 - 


10368X + 20736 


167 


( — .T^ 


— x^ — X^ — 1, X^ + X + 1) 


x^ 


- 4x5 


+ 2x4 _ 2x3 _ 3^,2 ^ 


h2x- 3 


177 


{x^ + 


x'^ + x^,x^ + x^ + 1) 




x^ 


+ 6x5 


+ 5x4 ^ g^3 _^ 2a;2 H 


hi 


188 




x'^ + x3 +x2 - 2x + 1,0) 




x5 


~xU 


h x3 + x2 - 2x + 1 




191 


( — X"^ 


+ x^ + x,x3 + X + 1) 




x^ 


+ 2x4 


- 2x3 ^ 5^2 g^, ^ 


1 



Table 4.1. Levels, integral models + h{x)y = ^(x), simplified models = /(x) 



Let us recall what is known about computing the quantities appearing on the right side of Equa- 
tion [T2] which we have not addressed so far. As described in [T3], the order of the torsion subgroups 
and the Tamagawa numbers are computable. For the Jacobians of the curves in Table 14. 1[ we 
list these values, taken from }131 Table 2], in Table |4?2] While no general algorithm has yet been 
developed and implemented to compute the order of the Shafarevich-Tate group in(A/(Q) for the 
Jacobians of each of the curves in Table 14.11 the conjectural order LU? of the group is also given, 
conditional on the classical BSD conjccture ll.il fand equal to the order of ni(A/Q)[2]). 

Remark 4.3. There is a general approach to computing L[I(A/Q) for the rank 2 Jacobians in Ta- 
ble which is to use Hcegner points and Kolyvagin's Euler system to give an explicit upper bound, 
then compute the remaining Selmer groups. It would be an interesting project to systematically de- 
velop this approach, by generalizing [121 123 to this new setting. 

Table l473l below provides the local height pairings away from p for iV 7^ 167. The global generators 
for A((Q)/A(Q)tors that we used are given as divisor classes [P-Q], [R-S], where P, Q, R, S e C{Q). 
Points at infinity are denoted by ooa, where a is equal to y/x^ evaluated at ooa- The heights list has 
three entries giving the nontrivial local height pairings hy{{P) — {Q), (i?) — (5)), hy{{P) — {Q), {—Q)~ 
(— P)) and hy{{R) — (S), {—S) — {—R)) for v ^ p. For two divisors D and E, this data is returned 
as a list of pairs [w, c?„], where v is a prime and hy{D, E) = dy ■ logp(w). 

Remark 4.4. The generators given for N ~ 125, A are actually generators for an index 2 subgroup 
of A(Q)/A(Q)tors, since an actual set of generators for the full group A{'Q)/A{Q)tors whose support 
solely consisted of non-Weierstrass points was not readily available. For N — 167, we had to use 
generators (of finite index subgroups) represented by divisors with pointwise Qp-rational support; 
see Section 1331 
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N 


Cy S 


|^((Q)tors| 


m? 


67 


1 


1 


1 


73 


1 


1 


1 


85 


4,2 


2 


1 


93 


4,1 


1 


1 


103 


1 


1 


1 


107 


1 


1 


1 


115 


4,1 


1 


1 


125,A 


1 


1 


1 


133,B 


1,1 


1 


1 


147 


2,2 


2 


1 


161 


4,1 


1 


1 


165 


4,2,2 


4 


1 


167 


1 


1 


1 


177 


1,1 


1 


1 


188 


9,1 


1 


1 


191 


1 


1 


1 



Table 4.2. BSD data for rank 2 Jacobians of genus 2 curves 



N 


global generators for A(Q)/A(Q)tors 


heights [[v, hy]] 


67 


[(-l,0)-oo_i],[(0,-l)-ooo] 


[],[],[] 


73 


[(-l,-2)-oo_i],[(0,-l)-ooo] 


[],[[3,1]],[] 


85 


[(-l,-2)-^_i],[(l,-4)-ooo] 


[[2,-l]],[[5,i]],[[5,i]] 
[]J[3,i]],[[3,i]] 


93 


[(-l,-2)-^_i],[(l,-3)-ooo] 


103 


[(-l,-l)-^_i],[(0,-l)-ooo] 


[],[],[] 


107 


[(-l,-l)-oo_i],[(l,-2)-ooo] 


[],[],[] 


115 


[(1,-4) [(-2,2) -ooo] 


[[3,-l]],[[5,i]],[[5,i]] 


125,A 


[(-l,0)-oo_i],[(l,-4)-ooo] 


[[2,-l]],[],[[5,l]] 


133,B 


[(0,-l)-oo_i],[(l,-2)-cx)o] 


[],[],[] 


147 


[(-l,-l)-^_i],[(-3,7)-ooo] 


[[2,-l]],[[3,i]],[[7,i]] 
[],[[7,i]],[[5,l],[7,|]] 


161 


[(l,-5)-oo_i],[(|,-3)-cx)o] 


165 


[(-8, -528) - (0,-144)], [(8, 80) - (0, 144)] 


[[2,2],[3,-i]], [[2,-2],[ll,i],[3,|]], 
[[2, -2], [5,1], [3, i]] 






177 


[(0,0)-oo_i],[(-f,-^)-cx)o] 


[[3,1]],[],[[3,-2],[17,1]] 


188 


[(0,l)-oo„i],[(-l,-l)-(2,5)] 


[[2,1]],[[2,|]],[[2,|],[5,1]] 


191 


[(0,-l)-oo_i],[(-2,10)-ooo] 


[],[], [[11,1]] 



Table 4.3. Global generators and intersectfon data 



For the computation of the special values £p{A,Q), we need the normalization factor 5~^, so 
we have to find a fundamental discriminant D > such that for some good ordinary pg we have 
gcd(;3o-^, D) = 1 and the analytic rank of is zero, where ip is the quadratic character associated 
to Q{VD). See Subsection El The real period for our Jacobians can be found in |13j : there it 
was computed using the observation that {uji,uj2) as in Remark l2.11l is a basis of integral 1-forms for 
all abelian varieties we consider. It is not difficult to show that the corresponding fact also holds for 
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all A^, and hence, using Remark 12.111 we found that 77^ = 1 in all 16 cases. We list, for each level, 
the quantities needed to find 5+ in Table 1131 



N 




'lip 


PO 








r (A 1 Vs 


I^VW/ torsi 


TIT / ? 


67 


5 


1 


19 


16 


4 


1/4 


1,1 


1 


4 


73 


5 


1 


11 


16 


4 


1/4 


1,1 


1 


4 


85 


61 


1 


41 


64 


16 


1/4 


4,2,8 


2 


1 


no 


r 



i 


ii 


— io 


4 


1 /A 

-1/4 


/111 

4,1,1 


i 


-1 
i 


103 


5 


1 


11 


16 


4 


1/4 


1,1 


1 


4 


107 


5 


1 


19 


16 


4 


1/4 


1,1 


1 


4 


115 


89 


1 


11 


64 


16 


1/4 


4,1,4 


1 


1 


125,A 


17 


1 


19 


16 


4 


1/4 


1,1 


1 


4 


133,B 


5 


1 


29 


-16 


4 


-1/4 


1,1,1 


1 


4 


147 


5 


1 


31 


-16 


4 


-1/4 


2,2,2 


2 


2 


161 


53 


1 


11 


64 


16 


1/4 


4,1,4 


1 


1 


165 


89 


1 


17 


64 


16 


1/4 


4,2,2,16 


4 


1 


167 


5 


1 


31 


-16 


4 


-1/4 


1,1 


1 


4 


177 


5 


1 


19 


16 


4 


1/4 


1,1,1 


1 


4 


188 


233 


1 


19 


144 


36 


1/4 


9,1,4 


1 


1 


191 


33 


1 


31 


16 


4 


1/4 


1,1 


1 


4 



Table 4.4. Rank zero twist data 



For good measure, we also verified: 

Proposition 4.5. The classical (and hence, for all primes p of good ordinary reduction, the p-adic) 
Birch and Swinnerton-Dyer conjecture holds for all 16 twists in Table under the assumption 
that III(A^/Q) is 2-torsion. 

Proof. Sec Table EH noting that m^,? is equal to |m(A^/Q)[2]|. □ 

Remark 4.6. For the computation of the Tamagawa numbers we used [71 Theorem 1.17]. Sup- 
pose that V I D, but that v does not divide the conductor of A. Then the twisted curve has 
bad reduction at v, but acquires good reduction over a quadratic extension. The classification of 
Namikawa and Ueno [28] shows that in this case must have reduction type [/q_o_o] ^-l ^- Since 
the geometric component group of the Ncron model is isomorphic to (Z/2Z)'^ for this reduction type, 
we always have c^{A^) \ 16. 

The tables below show the specific primes p and precision 0{p") for each level N for which we 
have tested Conjecture 11.41 

Remark 4.7. A note on our models and choices of primes. Since our p-adic heights algorithm 
requires that the curve be given by an odd degree model, for each curve = g{x), we consider those 
good ordinary primes p for which g(x) has a Qp-rational zero and do another change of coordinates 
to obtain the odd model y^ = f{x), with f{x) € Qp[a;]. We compute the p-adic regulators and p-adic 
L- values for these primes. 
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4.1. N ^ 67. We have the foUowing data: 



p-adic regulator Regp(A) 


p-adic L-value 


p-adic multiplier ep(A) 


905422 + 0(7«) 
655636176+ 0(13") 
6411910349+0(17**) 
1955457580+0(19^) 
6490501114813+ 0(37^) 
119112862323467+ 0(41^) 
231768637543452 + 0(43^) 
258343847102710 + 0(47^) 
7291679100956850+ 0(59^) 
6048812062982476+0(61^) 
53277934412195075+ 0(73^) 
9278983589215557+ 0(79^) 
157708559779041510+ 0(83^) 


4616447 + 0(7») 

3718847 + 0(13*) 

490126740 + 0(17*) 

205789013 + 0(19*) 

1520740814200 + 0(37*) 

6899026979535+0(41*) 

11662319738050+0(43*) 

6617527122585+0(47*) 

72703739307529+0(59*) 

174305066216353+0(61*) 

552479201354189+0(73*) 

88027589402068+0(79*) 

1578704504708054+0(83*) 


953283 + 0(7*) 

121846702 + 0(13") 

2996208382+0(17*) 

6722090086+0(19*) 

1763856795912+0(37*) 

604530321123+0(41*) 

11664993765232+0(43*) 

21577206386081 + 0(47*) 

77184936742982+0(59*) 

133406272889885+0(61*) 

460739420635942 + 0(73*) 

801037408797804+0(79*) 

162512920516158+0(83*) 



4.2. N = 73. We have the following data: 



p-adic regulator Regp(A) 


p-adic i-value 


p-adic multiplier ep(A) 


163731997+0(11*) 
482988818 + 0(13*) 
51174691892 + 0(23*) 
553299007790+0(31*) 
5421948177967+0(41*) 
38176784853304+0(59*) 
70602302343232+0(61*) 
433639741922965 + 0(71*) 
589304115938460+0(83*) 
6769596692483671 + 0(97*) 


183868925 + 0(11*) 
522787644+0(13*) 
46581832325 + 0(23*) 
440555494391 + 0(31*) 
5077531013725+0(41*) 
63020796753579+0(59*) 
139895606364222 + 0(61*) 
576931954734067+0(71*) 
347866087087015 + 0(83*) 
4269348271 + 0(97^) 


192773925 + 0(11») 

757562196 + 0(13*) 

48224542827+0(23*) 

850258335981 + 0(31*) 

5384419950679+0(41*) 

113039802800992+0(59*) 

79733480381568+0(61*) 

8989266238661 + 0(71*) 

720001059253854+0(83*) 

3993521258998096+0(97*) 



4.3. N = 85. We have the following data: 



j3-adic regulator Regp(A) 


p-adic L-value 


p-adic multiplier ep(^) 


1015073423894+0(37*) 

6819810980339+0(41*) 

8962714100713+0(53*) 

43568329449+0(61^) 

119416997911215+0(73*) 

1942338381733272 + 0(89*) 

5147606270477176+0(97*) 


167411116045+0(37*) 

7975900636623+0(41*) 

1069648287223+0(53^) 

3136884016567+0(61^) 

683019204724944+0(73*) 

3482744225118281 + 0(89*) 

2836855197 +0(97^) 


1002150510104+0(37*) 

7153783865856+0(41*) 

44683460285079+0(53*) 

72019680061615+0(61*) 

32602153132641 + 0(73*) 

1038134293650945+0(89*) 

3784121167774074+ 0(97*) 
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p-adic regulator Regp(A) 


p-adic L-vahic 


p-adic multipHer ep{A) 


185741420 + 0(11**) 
405483221 + 0(13^) 
43885519955 + 0(23®) 
336788503314+0(29**) 
1678741468628+0(37*) 
5324074002210+0(43**) 
16824305598488+0(47**) 
4960862919215 +0(538) 
143070222270789+0(61**) 
279322082363042 + 0(67**) 
430169136747961 + 0(738) 
1384453915387035 + 0(79^) 
2228621109604011 + 0(83^) 
1081659147745931 + 0(89**) 


151839057+0(11**) 
670790176+0(138) 
26161319539 + 0(23®) 
484038257980+0(29*) 
2569674002391 + 0(378) 
8725984878581 + 0(438) 
15669575995471 + 0(478) 
49317038818954 + 0(538) 
92506965732666+0(618) 
74111413499770 + 0(678) 
471513912864315+0(738) 
1232301086171477+0(798) 
1304530016876211 + 0(83*) 
1330994142123689+0(898) 


11051656 + 0(118) 
230057694 + 0(138) 
47769373949 + 0(23**) 
222828561623+ 0(298) 
1378603735422 + 0(378) 
3160248767946 + 0(438) 
8286455636222 + 0(478) 
23178143892193 + 0(538) 
122923764909639+0(618) 
180555937022120+ 0(678) 
411163460552347+0(738) 
503819572894975+ 0(798) 
1071794429225898+ 0(838) 
1518075886594725+ 0(898) 



4.5. N = 103. We have the followmg data: 



p-adic regulator Regp(A) 


p-adic i-value 


p-adic multiplier f-p{A) 


147377758 + 0(118) 
489193484 + 0(138) 
15204606664+0(198) 
66216995216 + 0(238) 
5372718408 + 0(298) 
6799091682040+0(418) 
23467041445332 + 0(478) 
9449958206985+0(538) 
55788681659810+ 0(598) 
180708198470076+0(618) 
304798054862709+0(718) 
651632900917334+0(738) 
1422073004111088+0(798) 
2204776989584744+0(838) 
3419478873681093+0(898) 
5187359554281130 + 0(978) 


86486502 + 0(118) 
67428377 + 0(138) 
10638300382 + 0(198) 
18109392006+0(238) 
76731347688+0(298) 
4391281006909+0(418) 
5937816898560 + 0(478) 
6585582284426 + 0(538) 
58416917952322+0(598) 
86273076603078+0(618) 
23644536785282+0(718) 
128186925484 +0(736) 
162819364440040 + 0(798) 
290525162365+0(83*^) 
13722104837501 + 0(89^) 
588713923936+0(97*^) 


192773925 + 0(118) 

337691501 + 0(138) 

12173049603+0(198) 

45043095109+0(238) 

12536647436+0(298) 

5303120857798+0(418) 

1847891549858 + 0(478) 

47071170371848+0(538) 

134523844728309+0(598) 

91320952633362+0(618) 

318837731560077+0(718) 

86431680403618 + 0(738) 

1353067258168647+0(798) 

1489257165084816+0(838) 

2011257469143583 + 0(898) 

291254315420391 + 0(978) 
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4.6. 



N = 
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107. We have the fohowing data: 



p-adic regulator Regp(A) 


p-adic L- value 


p-adic multiplier £p{A) 


100037184 +0(13«) 
5164824485 + 0(17^) 
4948122310 + 0(19^) 
2233155353996 +0(37«) 
7693933727093+0(418) 
1985728518871 + 0(438) 
9236325503676 +0(478) 
66592510503713+0(598) 
73605475872145+0(618) 
215631855830774+ 0(678) 
235988007934369+0(718) 
1405009654786451 + 0(798) 
966246807067004+0(838) 


381034778 + 0(138) 
2251756830+0(178) 
410682533 + 0(198) 
782254360600 + 0(378) 
1126806110759+0(418) 
11350309489829+0(438) 
9352977397857 + 0(478) 
126498662012390+0(598) 
67998854641813+0(618) 
276136144242399 + 0(678) 
4203263257242+0(71^) 
1425538781612665+ 0(798) 
29073683 +0(83^) 


725504508+0(138) 
6548185060 + 0(178) 
11770828305+0(198) 
499167048517 + 0(378) 
941047246375+0(418) 
1490080422844 + 0(438) 
3877324130340 + 0(478) 
66521146158463+0(598) 
149495311709314+0(618) 
300204582979356 + 0(678) 
297538516047502+0(718) 
637219753066297+0(798) 
122301722091732 + 0(838) 


N = 115. We have the following data: 


p-adic regulator Regp(A) 


p-adic L-value 


p-adic multiplier ep{A) 


151819184+0(118) 
6070540659 + 0(178) 
443043366998 + 0(378) 
10506890337861 + 0(438) 
25938666299194+ 0(538) 
68828469915327+0(598) 
117125015025879+0(618) 
117261157211649 + 0(678) 
38346420175144+0(798) 
232154244720909+0(838) 
3680613169329886+0(898) 
337111037730418 + 0(978) 


96031694+ 0(11») 
2602174031 + 0(178) 
3207000318071 + 0(378) 
466034248434+0(438) 
43404652273198+0(538) 
15822514736163+0(598) 
99513360280408+0(618) 
369365142758789 + 0(678) 
1042744621946608+0(798) 
1238796074898239 + 0(838) 
1982664616252635 + 0(898) 
5523549952859660 + 0(978) 


55470083 + 0(118) 
2479111430+0(178) 
518402902203 + 0(378) 
3160248767946 + 0(438) 
24704105954182 + 0(538) 
112268718282797+ 0(598) 
144234417021077+0(618) 
388573100762289+ 0(678) 
856522414733559+0(798) 
2232922964727286+ 0(838) 
1447736508567520 + 0(898) 
5537610452725212+0(978) 


N — 125, A. We have the following data: 


p-adic regulator Regp(A) 


p-adic i- value 


p-adic multiplier ep{A) 


298562498 + 0(138) 
6712555657+0(198) 
28761182485+0(238) 
1610334394992 + 0(378) 
5827125727855 + 0(478) 
43827481404730 + 0(538) 
934069839446 + 0(598) 
94940897306587+0(618) 
344652595573416 + 0(678) 
494778091759992 + 0(738) 
911058348384486+0(838) 
3812663593637783+0(898) 
7743507247513256 + 0(978) 


592894408 + 0(138) 

7153379737+0(198) 

19244567041 + 0(238) 

2619837199442 + 0(378) 

2438975823319 + 0(478) 

12237145144494+0(538) 

138331812786050+0(598) 

86820693223899+0(618) 

285200220171958 + 0(678) 

211635131306627+0(738) 

15062127863580+0(83^) 

140584030153+0(89*^) 

72751227747664+0(97^) 


337691501 + 0(138) 
7352726322 + 0(198) 
47769373949+0(238) 
518402902203 + 0(378) 
20137488978024 + 0(478) 
45444073360562 + 0(538) 
119463235911829+0(598) 
134545325721836+0(618) 
55395450190703 + 0(678) 
338214791799846+0(738) 
1695835531921770+0(838) 
3229221353736449+0(898) 
6602401135477806 + 0(978) 



4.9. N = 
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133,5. We have the foUowmg data: 
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p-adic regulator Regp(A) 


p-adic L-value 


p-adic multiplier ep{A) 


4554714851 + 0(17«) 
482641533 +0(296) 
285247284517+ 0(31^) 
873461875052+0(418) 
6395433286380+0(43*^) 
40174155934745+0(53*) 
388303423009987+ 0(67*) 
582542046575002 + 0(73*) 
997934987934019+0(79*) 
337083794306147+ 0(83*) 
3826161118964265+0(89*) 


1400369830 + 0(17*) 
224834369110+0(29*) 
644745508559+0(31*) 
5913841764921 + 0(41*) 
5250591893580+0(43*) 
9436443664 +0(536) 
74275805470 +0(676) 
47404160292+0(736) 
1692929332309+0(79^) 
8815903470 + 0(836) 
689438763 + 0(89^) 


1530767973 + 0(17*) 

188246220652+0(29*) 

65426082523+0(31*) 

5173622706020+0(41*) 

7173815953060+0(43*) 

47071170371848+0(53*) 

287829738202699+0(67*) 

214334244118640+0(73*) 

85649658584845+0(79*) 

1446666792043837+0(83*) 

1571471061650586+0(89*) 


N = 147. We have the following data: 


p-adic regulator Regp(A) 


p-adic i- value 


p-adic multiplier ep{A) 


434194800 + 0(13*) 

3085885399 + 0(19*) 

57105870 + 0(23^) 

598807495296 + 0(31*) 

1255556858069+0(37*) 

3028914438423+0(43*) 

21722415097178+0(53*) 

9696531871680+0(61*) 

252460432397529+0(67*) 

319985315705867+0(71*) 

696485497462517+0(73*) 

1036811888178773+0(79*) 

4273472549945572 + 0(97*) 


772553365 + 0(13*) 
14351355419+0(19*) 
56314647135 + 0(23*) 
500890389807+0(31*) 
728534804896+0(37*) 
5310811645878+0(43*) 
32574036544128 + 0(53*) 
147051772023912 + 0(61*) 
2933215623449+0(67^) 
4744056079140+0(71^) 
75645384726+0(736) 
965207536 +0(79^) 
6770845150 +0(97^) 


69777210 + 0(13*) 

10124513344+0(19*) 

77688619426+0(23*) 

226083261470+0(31*) 

925099803678+0(37*) 

1448245155768 + 0(43*) 

14098536063957+0(53*) 

127085340697404+0(61*) 

15805729099128+0(67*) 

454718387048106+0(71*) 

263081220212640+0(73*) 

74020902743243+0(79*) 

7422648274246094+0(97*) 


N — 161. We have the following data: 


p-adic regulator Regp(A) 


p-adic L-value 


p-adic multiplier ep{A) 


171933135 + 0(11*) 
16676191757+0(19*) 
3000539180980+0(37*) 
4799012913812+0(43*) 
33038825747471 + 0(53*) 
144048375212404+0(59*) 
989091293021 + 0(61*) 
93625125465306+0(67*) 
372742847896101 + 0(79*) 
133689266642605 + 0(83*) 
109602346601919+0(89*) 
4449889265258731 + 0(97*) 


104178769 + 0(11*) 
8396822512+0(19*) 
2959738471101 + 0(37*) 
820015420 + 0(436) 
14752347930+0(536) 
67829338510607+0(59*) 
62763431617869+0(61*) 
43269028077 + 0(676) 
16260012523515+0(79^) 
27328857470 + 0(836) 
3775670578 +0(89^) 
8056056109 +0(97^) 


48803991 + 0(11*) 

10186228540+0(19*) 

1378603735422+0(37*) 

7358928540810+0(43*) 

37415160388754+0(53*) 

113039802800992+0(59*) 

191252449121304+0(61*) 

228681540167106+0(67*) 

1380060506871347+0(79*) 

1997163923487638 + 0(83*) 

404117712562583+0(89*) 

3796862465389610 + 0(97*) 
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165. We have the fohowing data: 



p-adic regulator Regp(A) 


p-adic L- value 


p-adic multiplier ep{A) 


2478665 + 0(7^^) 
7577669996 + 0(13^) 
1345832 + 0(17^) 
317314039860+0(19^) 
1197529401836+0(23^) 
182820405709 +0(29^) 
30402585606264+ 0(37^) 
300867423531184+ 0(41^) 
453491841293220+ 0(43^) 
841500704008115 +0(47^) 
108480654690546+ 0(53^) 
5924946879989069 + 0(59^) 
11004676059690151 + 0(61^) 
6402333518135195 + 0(67^) 
18266464992713450+ 0(71^) 
42183534718264644+ 0(73^) 
1224455456912234+ 0(798) 
124036966428761339+ 0(83^) 
243696118400513337+ 0(89^) 
635540819872824429+ 0(97^) 


988615 + 0(7») 
546478360+0(13*^) 
115518752 + 0(17^) 
15454527827+ 0(198) 
21430827992 + 0(238) 
462793840863+0(298) 
145763317789 + 0(378) 
4638175450295+0(418) 
10502890759714 + 0(438) 
22177682954670 + 0(478) 
812837848921 + 0(53^) 
2021593887077+0(59^) 
1114772875983+0(61^) 
1539009404714+0(67^) 
98457973781 + 0(71^) 
6514036760733+0(73^) 
1812891052 + 0(79^) 
203487283131 + 0(836) 
2097535192 + 0(89^) 
1063985237 +0(97^) 


2047938 + 0(78) 
19120487+0(138) 
6743866153 + 0(178) 
16701261693 + 0(198) 
32283075894 + 0(238) 
167834932128+0(298) 
1905855970461 + 0(378) 
7243944162192 + 0(418) 
981850330755 + 0(438) 
15840901508219 + 0(478) 
17925543534180 + 0(538) 
32899436516884+0(598) 
102536167075224+ 0(618) 
15805729099128 + 0(678) 
34536869719889+0(718) 
39739465931437 + 0(738) 
186442021878008 + 0(798) 
1504086404024377 + 0(838) 
1826662988317474 + 0(898) 
4391966065852909+ 0(978) 



4.13. N = 167. We have the following data: 



p-adic regulator Regp(A) 


p-adic L- value 


p-adic multiplier ep{A) 


19432714 + 0(7'-') 
13117611 + 0(13^) 
1908862518313+0(171°) 
2053675284265 + 0(1910) 
27719111127295 +0(231°) 
405640880151858+0(311°) 
57769565310991429 + 0(53i°) 
501328316424338015+0(591") 
1243828341260907954+ 0(7li°) 
2079988387733147685+0(731°) 
20365783254113182401 + 0(89i°) 
67990777180272953115+0(971°) 


2251 + 0(7") 
666390377 + 0(138) 
2314174880 + 0(178) 
7656154501 + 0(198) 
41736439730 + 0(238) 
714765172682 + 0(318) 
19704952386 +0(536) 
64487204069600+0(598) 
56573288611 + 0(7l6) 
4104591 + 0(73"^) 
5065696436 + 0(89^) 
23364634 + 0(97'') 


307185 + 0(78) 
526042526+0(138) 
5234654956+0(178) 
14340680958 + 0(198) 
16647712571 + 0(238) 
358709025654+0(318) 
17681237786119+0(538) 
140344643451642 + 0(598) 
557648531014830+0(718) 
759308640111719+0(738) 
1447017073110591 + 0(898) 
6376229493766338 + 0(978) 
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4.14. N = 177. We have the fohowing data: 
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p-adic regulator ticg^(A) 


p-adic _L-value 


p-adic multiplier ep[A) 


1072267 +0(7») 
9772408 + 0(196) 
27690468499+0(23®) 
141718660962+0(29®) 
265288097732 + 0(31®) 
1019326123826+0(37®) 
952644023485 + 0(41®) 
21867793727731 + 0(47®) 
54813744728211 + 0(61®) 
374976464608823+0(73®) 
2024750045809193+0(83®) 


1192 + 0(7*) 
2558009183+0(19®) 
51308343838+0(23®) 
391909937451 + 0(29®) 
167635394515 + 0(31®) 
2370016933013+0(37®) 
5888249521909+0(41®) 
3399186192198+0(47®) 
136070004398022 + 0(61®) 
92506712920 +0(736) 
1537698302 + 0(836) 


507488 + 0(7®) 

11268267357+0(19®) 

3788873485 + 0(23®) 

65127401733+0(29®) 

519021947371 + 0(31®) 

1021993916814+0(37®) 

5919384948361 + 0(41®) 

18039627327365+0(47®) 

149495311709314+0(61®) 

469524064138469+0(73®) 

2034477952337988 + 0(83®) 



4.15. N = 188. We have the following data: 



p-adic regulator Regp(A) 


p-adic L-value 


p-adie multiplier ep{A) 


5623044+0(7®) 

4478725 + 0(11^) 

775568547+0(13®) 

1129909080+0(17®) 

14409374565+0(19®) 

31414366115+0(23®) 

2114154456754+0(37®) 

6279643012659+0(41®) 

9585122287133+0(43®) 

3328142761956 + 0(53®) 

17411023818285+0(59®) 

102563258757138+0(61®) 

26014679325501 + 0(67®) 

490864897182147+0(71®) 

689452389265311 + 0(73®) 

878760549863821 + 0(79®) 

2070648686579466+0(83®) 

3431343284115672 + 0(89®) 

4259144286293285 + 0(97®) 


1259 + 0(7*) 
150222285 + 0(11®) 
237088204+0(13®) 
6922098082 + 0(17®) 
15793371104+0(19®) 
210465118 + 0(23®) 
1652087821140+0(37®) 
2066767021277+0(41®) 
3309737400961 + 0(43®) 
5143002859 + 0(536) 
7961878705 + 0(596) 
216695090848 + 0(61^) 
7767410995 + 0(676) 
16754252742 + 0(716) 
193236387 +0(73^) 
1745712500 + 0(79^) 
2888081539 + 0(83^) 
1591745960 + 0(89^) 
21828881 + 0(97*) 


507488 + 0(7®) 

143254320 + 0(11®) 

523887415 + 0(13®) 

4494443586 + 0(17®) 

4742010391 + 0(19®) 

45043095109+0(23®) 

1881820314237+0(37®) 

4367414685819+0(41®) 

85925017348+0(43®) 

6112104707558+0(53®) 

98405729721193+0(59®) 

137187998566490+0(61®) 

38320151289262+0(67®) 

530974572239623+0(71®) 

162807895476311 + 0(73®) 

1063642669147985 + 0(79®) 

1103760059074178+0(83®) 

1012791564080640+0(89®) 

6376229493766338+0(97®) 
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4.16. N = 191. We have the fohowing data: 



p-adic regulator Rcgp(j4) 


p-adic L- value 


p-adic multiplier ep{A) 


4195478 + 0(7«) 
43495803539 +0(238) 
276478270993 +0(31*) 
7847912037610 + 0(43**) 
3701160666066+0(47*) 
19837992635361 + 0(53*) 
207820830309704+0(71*) 
105659818394179 + 0(73*) 
4330286071100495 + 0(97*) 


1867 + 0(7-*) 
62365909362 + 0(23*) 
411081898951 + 0(31*) 
1839263047933+0(43*) 
16594732090932+0(47*) 
121641372 + 0(535) 
80098460638243+0(71*) 
278456920 + 0(73^) 
12214648 + 0(97") 


1638463 + 0(7*) 

47598354917+0(23*) 

611200443823+0(31*) 

10085036614653+0(43*) 

9836262988784+0(47*) 

22289116823061 + 0(53*) 

318837731560077+0(71*) 

160255667550084+0(73*) 

1683523428082670 + 0(97*) 



5. Evidence for a twist of rank 4 

In this section, we present evidence for Conjecture 11.41 on a rank 4 twist of a rank modular 
abelian surface for the primes 29, 61 and 79. Let X = Xo{31). According to [T?], an affine equation 
for X is given by 

= {x^ - 2^2 - X + 3) • {x^ - 6x^ - 5x - 1). 
The Jacobian A oi X has rank zero over Q. 

We search for quadratic twists of high rank by searching among quadratic twists ip of small 
conductor D for some A-^ whose complex L-series seems to vanish at s = 1 up to order at least 4. 
This is the case for D = —47. We then use a 2-descent on A^p as implemented in Magma to find 
that the rank is at most 4. Searching for Q-rational points on A^ of small height quickly reveals 
subgroups of rank 4, such as the groups Gp described below, thus proving that the rank A^p over Q 
is indeed 4. 

Using [71 Theorem 1.17] we find that the potentially nontrivial Tamagawa numbers arc C31 {A^) = 1 
and C47(A^) = 16. Moreover, the torsion subgroup is trivial as is the 2-torsion of III(A^,/Q). Since 
the divisors supported in Q-rational points of X^ do not generate a subgroup of finite index, we 
compute the p-adic regulator for p G {29,61, 79} using the second method outlined in Section 13.41 
Namely, for each p we find a finite index subgroup Gp of A^(Q)/A^(Q)tors generated by the classes 
of divisors Di, . . . , Z?4 such that each Di has pointwise Qp-rational support. It is then enough to 
compute the p-adic regulator of Gp to find the p-adic regulator of A.ff,{<Q)/A,p{Q)toi-s- 

The generators we used are given in Table 15.11 in Mumford representation, along with the index 
of G„. 



p index generators of Gp 

[x-^ - 7/2a- + 49/16, 581/16a; - 1305/32], [x^ - 2x - 1/2, -47/2a;], 
[x^ - 5/3a; + 5/6, 47/18a; - 517/18], [x^ - 19/3a; - 11/4, -517/36x - 47/12] 
61 4 - 7/2a; + 49/16, 581/16X- 1305/32], [x^ - 5a; + 11/2, -235/2x + 423/2], 

[x^ - 5x- 7/3, 235/3a; + 94/3], [x^ - 5/3x + 5/6, -47/18a; + 517/18] 
79 4 [x2-7/2a; + 49/16,581/162;- 1305/32], 1/3, -47/3a;], 

I I [x^ - 3a; - 5/3, 149/3a; + 82/3], [x^ - 19/7x + 3/7, -1363/49a; - 564/49] 

Table 5.1. Indices and generators, = 31 twisted hy D = —47 

In order to compute the special values of the p-adic L-series, we need to find out the correct 
normalization factor of the minus modular symbol map associated to A. The twist A^ of A by 



p-ADIC BSD FOR MODULAR ABELIAN VARIETIES 31 



the quadratic character x associated to Q(\/~19) has rank over Q. Coniparing [0]j^^ = —4 to 
^^^^^ = 1, we find that $- = -i. 

For the latter computation, we used that for a minimal equation of X^, a basis of the integral 
1-forms on is given by {uji,uj2) as in Remark l2.11l Since the corresponding fact also holds for A, 
we find that % = 1 using [TO] . 

The data presented in this section proves: 

Proposition 5.1. Assume that III(A^/,/Q) consists entirely of 2-torsion, and that the conjectural 
order o/III(yl^/Q) is 1 (numerically it is 1.0000000... to as many digits as we care to compute). 
Then Conjecture \1.4\ is satisfied for the twist A^, of Jq(31) of rank 4 for the primes 29, 61 and 79. 

The special values of the p-adic L-series, the p-adic regulators and the p-adic multipliers for 
p € {29, 61, 79} are given in the following table: 



p-adic regulator Regp(A) 


p-adic L- value 


p-adic multiplier ep{A) 


351486231941615978+0(29^^) 

1650697608489237057465+0(6112) 

8155329946924028539010+0(79^2) 


202402009906 +0(29«) 
4326648666405 + 0(61**) 
1513185184992411+ 0(798) 


423952915488 +0(29») 
10267186717780+0(61*) 
1431106352547896 + 0(79^) 



Remark 5.2. Assume that III(74^/Q) consists entirely of 2-torsion. Then we also verified the 
classical (and hence for good ordinary primes p the p-adic) conjecture of Birch and Swinnerton-Dyer 
for the rank twist A^, since all Tamagawa numbers, the order of the torsion subgroup and the 
order of III(A;^/Q)[2] are easily seen to be equal to 1. 
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